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ABSTRACT 


A  comparison  of  a  high  speed  container  ship  using  a  rudder  versus  a  podded 
propulsor  is  made  to  study  replacing  a  rudder  with  a  pod.  A  mathematical  model  is 
altered  to  simulate  a  ship  operating  with  a  rudder  and  with  a  pod  to  maneuver.  The  model 
incorporates  the  nonlinear  maneuvering  equations  and  couples  the  surge  and  sway  forces, 
yaw  and  roll  moment,  and  the  roll  angle  induced  during  a  steady  turn  with  varying  rudder 
and  pod  angles.  The  model  uses  the  hydrodynamic  derivatives  and  coefficients  for  a  high 
speed  container  ship.  The  equations  are  numerically  integrated  in  order  to  predict  the  roll 
angle,  sway  and  surge  velocities,  and  the  ship’s  position  in  the  xy-plane.  Both  transient 
and  steady  state  results  are  utilized  to  quantify  the  relative  efficiency  of  each  system.  The 
results  are  used  as  a  preliminary  study  into  replacing  a  rudder  on  a  ship  with  a  podded 
propulsor.  The  results  indicate  that  the  ship  responds  faster  and  has  a  shorter  turning 
radius  with  the  pod  at  lower  initial  speeds  and  pod  angles,  while  the  rudder  responds 
better  at  high  speeds  regardless  of  angle.  Further  research  is  necessary  to  study  the  effects 
of  changing  the  pod’s  position  and  increasing  the  number  of  pods  used. 
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I.  INTRODUCTION 


A.  GENERAL 

Most  modern  ships  use  a  rudder  to  provide  the  necessary  forces  and  moments  to 
allow  the  ship  to  turn.  However,  in  recent  years  the  development  of  podded  propulsors 
has  allowed  ships  with  pods  to  increase  their  propulsive  efficiency  and  maneuverability 
[Ref  1].  The  use  of  pods  allows  the  ship  to  operate  with  a  diesel-electric  power  plant.  If 
the  pods  are  aft  on  the  ship,  the  pod  can  increase  the  forces  generated  while  maneuvering 
at  low  speeds  and  eliminate  the  need  for  stern  thrusters  [Ref.  1].  Podded  propulsors  are 
currently  used  in  cruise  ships,  where  passenger  comfort  is  the  main  concern,  dredging 
vessels,  and  naval  vessels  such  as  mine  hunters.  The  pod  allows  the  ship  to  maneuver 
more  rapidly  with  less  adverse  effect  on  the  ship  during  turns,  especially  at  slow  speeds. 
Thus  the  ship  can  turn  faster  at  low  speeds  and  will  have  smaller  roll  angles.  The  ship 
also  has  a  smaller  turning  diameter  at  slow  speeds  and  pod  deflections  [Ref.  2].  Another 
advantage  of  pods  is  that  they  can  be  steered  over  360°,  allowing  the  ship  to  maneuver  in 
any  direction  fairly  easily.  Pods  also  take  up  less  space  on  a  ship  because  the  need  for 
long  shaft-lines,  conventional  rudders,  and  reduction  gears  is  eliminated  [Ref.  2]. 

Pods  operate  in  a  similar  fashion  as  rudders.  To  turn  a  ship,  one  simply  positions 
the  rudder  at  the  desired  angle  to  accomplish  the  turn  at  a  certain  rate.  A  ship  with  a  pod 
can  apply  the  necessary  angle  by  rotating  one  or  all  of  the  pods,  if  more  than  one  pod  is 
installed  [Ref.  1].  Because  of  the  weight  of  the  pod,  turning  by  using  a  pod  is  considered 
to  cause  more  strain  on  the  equipment  than  steering  with  rudders.  However,  rudders 
experience  stall  at  about  35  degrees,  whereas  pods  do  not  experience  stall  at  all  due  to  the 
fact  that  the  force  is  in  the  same  direction  as  the  propeller  thrust  [Ref.  1].  Rudders  are 
also  used  to  help  maintain  the  course  of  a  ship  when  attempting  to  travel  in  a  straight  line. 
The  ship  may  veer  from  its  intended  course  due  to  ocean  waves  and  currents.  According 
to  Toxopeus  and  Loeff,  vessels  with  podded  propulsion  perform  well  at  course  keeping 
and  yaw  checking  [Ref.  1]. 

However,  pods  do  have  some  disadvantages.  Pods  can  respond  very  rapidly  and 
therefore  may  respond  faster  than  a  human  operator,  thus  making  it  difficult  for  a 
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helmsman  to  control  the  ship  [Ref.  2],  One  solution  may  be  to  implement  computer 
controls  to  allocate  the  thrusts  of  the  propeller  in  such  a  way  that  the  environmental  loads 
do  not  have  a  detrimental  effect  on  the  ship’s  perfonnance.  Pods  also  tend  to  cause  large 
heel  angles.  According  to  a  study  performed  by  MARIN,  most  passengers  tend  to  panic  at 
heel  angles  of  7°,  and  pods  have  exhibited  heel  angles  of  up  to  25°  in  some  ship  models 
[Ref.  2].  Also,  in  order  for  the  pod  to  operate  in  the  most  effective  manner  possible,  the 
ship  hull  has  to  be  altered  and  usually  more  than  one  pod  is  required  to  provide  adequate 
power  and  maneuverability  [Ref.  1].  However,  because  pods  have  proven  so  versatile, 
further  study  should  be  done  to  see  if  they  can  be  used  in  more  applications. 

B.  SCOPE  OF  THIS  THESIS 

The  purpose  of  this  thesis  is  to  compare  the  performance  of  a  ship  using  a  pod 
with  that  of  a  ship  using  a  rudder  with  different  speeds  and  rudder  and  pod  deflections. 
The  comparison  is  made  to  help  determine  the  viability  of  replacing  rudders  on  naval 
vessels  with  pods.  Since  naval  vessels  may  need  to  maneuver  rapidly  at  high  speeds,  for 
example  to  avoid  a  torpedo,  it  is  important  to  detennine  if  a  pod  will  perform  well  under 
such  conditions. 

The  comparison  of  the  two  ships  is  made  using  a  model  that  incorporates  data  for 
a  high  speed  container  ship.  The  model,  described  in  Chapter  II,  uses  the  equations  of 
motion  for  a  ship  to  develop  equations  for  the  derivatives  of  several  parameters.  The 
original  model  simulates  a  ship  with  a  rudder.  The  model  is  modified  to  simulate  the  use 
of  a  podded  propulsor  rather  than  a  rudder.  It  is  important  to  note  that  the  modifications 
made  on  the  model  do  not  take  into  account  any  change  in  the  hull  that  may  be  required 
for  use  of  a  pod,  nor  the  possible  need  of  more  than  one  pod  to  properly  maneuver.  The 
model  simply  removes  the  rudder  and  assumes  that  the  propeller,  whose  characteristics 
have  not  been  altered,  can  now  rotate  in  the  same  manner  as  a  pod. 

The  analysis  performed  will  help  determine  the  general  operating  characteristics 
of  a  podded  propulsor  over  a  range  of  speeds  and  pod  deflections.  A  description  of  the 
transient  analysis  perfonned  is  found  in  Chapter  III.  The  steady  state  analysis  performed 
is  described  in  Chapter  IV.  The  parameters  for  each  ship  that  are  compared  are:  the  total 
speed  U,  the  sway  velocity  v,  the  roll  angle  (j),  and  the  ship’s  position  in  the  xy-plane. 
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These  parameters  are  plotted  for  each  ship  configuration  and  are  then  used  to  study  the 
pod’s  operating  characteristics. 

The  conclusions  drawn  from  the  analysis  performed  are  shown  in  Chapter  V. 
Areas  of  further  research  are  also  recommended. 
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II.  MATHEMATICAL  MODEL 


A.  EQUATIONS  OF  MOTION 


The  nonlinear  equations  of  motion  in  a  body-fixed  coordinate  system  for  a  ship 
are  shown  in  equations  (2.1).  The  equations  assume  symmetry  along  the  longitudinal  axis 
(yG=0).  The  orthogonal,  right-hand  coordinate  system  is  shown  in  Figure  1  and  the 
parameters  are  defined  in  Table  I. 

Surge:  myu -xG  (q2  +r2 )  +  zG  (q  +  pr}-rv  +  wqj  =  X  (2.1. a) 

Sway:  m f v  +  xG  ( pq  +  r )  +  zG(qr  -  p)  +  ur-wp\  =  Y  (2.1.b) 


Surge:  m^u  -xG(q2  +r2  )  +  zG  (q  +  pr)-  rv  +  wqj  =  X 

Sway:  m  [  v  +  xG  ( pq  +  r  )  +  zG  ( qr  -  p  )  +  ur  -  wp ]  =  Y 

Heave:  m  ^w  +  xG  ( rp  -  q)- zG(p2  +  q2^j  +  pv  -qu^  =  Z 

n  „  IxP  +  {l;-I,)qr-Ixy{pq  +  f)-mzG{v  +  ur-wp)  = 
Roll: 

K-AGZ(<f>) 

D.  ,  IyA-(l:-Ixyp-Ix\r2-p2)  + 

Pitch: 

myzG  (u  +  qw-i~v^-xG  (w  +  pv-  qu)~^  =  M 
Yaw:  Ij  +  (ly  - 1, )  pq  +  Ixz  (qr  -  p)  +  mxG  (v  +  ur  -wp)  =  N 


(2.1.c) 


(2-1-d) 


(2.1.e) 


(2-l.f) 


Figure  1.  Coordinate  System  [Ref.  7] 
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Table  1  Equations  of  Motion:  Parameters 


Parameter 

Description 

x,  y,  z 

Distance  along  the  principal  axes 

U,  V,  w 

Translational  velocity  components  of  ship  relative  to  fluid 
along  body  axes 

P,  q.r 

Rotational  velocity  components  of  ship  relative  to  inertial 
reference  system  along  body  axes 

i  X,  Y,  Z 

Hydrodynamic  force  components  along  body  axes 

K,  M,  N 

Hydrodynamic  moment  components  along  body  axes 

M V,  0,  (|> 

v| /  yaw  angle:  bow  to  starboard  positive 

0  pitch  angle:  bow  up  positive 

cj)  roll  angle:  starboard  down  positive 

m 

Mass  of  ship 

xg,  Yg,  zg 

Coordinates  of  the  center  of  gravity  in  the  body  axis 
system 

lx,  Iy,  Iz 

Moments  of  inertia  about  the  body  axis  system 

Ixz,  Iyz,  Ixy 

Products  of  inertia  about  the  body  axis  system 

V 

Displacement  volume  of  ship 

!  A 

Displacement  weight  of  ship 

1  GZ(d>) 

Righting  moment  as  a  function  of  roll  angle 

1  5 

Rudder  angle  in  radians 

V 

Initial  velocity  of  ship 

p 

Mass  density  of  sea  water 

L 

Ship  length  between  perpendiculars  (LBP) 

B.  SIMPLIFYING  ASSUMPTIONS 

To  simulate  a  ship  maneuvering  in  the  ocean  without  changes  in  depth,  the 
following  assumptions  were  made  [Ref.  3]: 

1.  The  rotational  velocity  and  acceleration  about  the  y-axis  are  zero,  (q  =  0 
and  q  =  0  ) 

2.  The  translational  velocity  and  acceleration  in  the  z  direction  are  zero. 
(w  =  0  and  w  =  0 ) 

3.  The  vertical  heave  and  pitch  motions  are  decoupled  from  the  horizontal 
plane  motions. 

4.  The  product  of  inertia  IXz  is  small  and  can  be  neglected. 
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5.  The  surge  equation  is  substituted  by  an  algebraic  equation  which  is  a 
function  of  u,  V,  and  5. 

6.  The  longitudinal  center  of  gravity,  (LCG),  and  the  longitudinal  center  of 
buoyancy,  (LCB),  are  at  midship. 

7.  The  vertical  center  of  gravity,  (VCG),  is  on  the  centerline. 

8.  The  only  important  forces  and  moments  acting  on  the  ship  induced  by  the 
rudder  (or  pod)  are  due  to  rudder  (or  pod)  deflection.  Forces  and  moments 
due  to  8  and  8  are  negligible.  [Ref.  4] 

Applying  these  assumptions  to  Equations  (2.1),  the  following  simplified  equations 
of  motion  are  obtained  [Ref.  5]: 


Sway:  m  [v  +  xGr  -  zGp  +  ur 1  =  Y 

(2. 2. a) 

Roll:  Ixp  -  mzG  (v  +  ur)  =  K  -  AGZ  ((f)) 

(2.2. b) 

Yaw:  I.r  +  mxG  ( v  +  ur )  =  N 

(2.2.c) 

To  simplify  the  analysis  even  further,  equations  (2.2)  are  converted  to  a 
nondimensional  form.  The  nondimensional  relationships  are  shown  in  Table  2. 

To  demonstrate  the  conversion  to  nondimensional  tenns,  the  appropriate  values 
from  Table  2  are  substituted  in  Equation  (2.2a)  to  yield  the  following: 


\pl* 

2  j 


in 


f  y2\ 
V  L  J 


v  +  Lx, 


V2 

T2 


-\Lzr 


IL 


+  (Vu) 


V  , 

—  r 
L 


f\  9  2 

-pL2V~ 


r 


Simplifying  by  factoring  out  (V“/L)  gives 


f  1 

-pllv2 

v2 


lit 


•  f  ,  f  • !  f  •  t  .  ft 

v  +xGr  -zG  p  +ur 


f  1 

-pL2V2 

v2 


r 


Equation  (2.2a)  can  then  be  rewritten  as  Equation  (2.3)  as  shown. 


m 


•  f  .  r  •  t  f  •  f  ,  ft 

v  +xGr  —zG  p  +ur 


=  r 


(2.3) 
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Because  Equations  (2.2)  and  (2.3)  are  of  the  same  form,  the  prime  notation  will  no  longer 
be  included  and  all  equations  will  be  treated  as  nondimensional  unless  otherwise 
specified. 

Table  2  Nondimensional  Parameter  Relationships 


,  V 

v  =  — 

V 

/ 

v'  =  v 

V 

L  > 

V2, 

Y'  =  Y 

f  1  ^ 

(1/2  Pi2vy 

t 

r  —r 

'  L  \ 

f 

•  t 

r  =  r 

V 

L2^ 

N'  =  N 

(  1 

\ 

/ 

U2J 

^1/2  plPV2 

ii 

II 

(L' 

K'  =  K 

1 

S/ipilv 2J 

p  =  p 

' L N 

P'  =  P 

f  L2 ' 

y2. 

i  u 

u  =  — 

V 

m  =  m 

<  1  ^ 

,1/2  pi5  J 

XG  =  *< 

(i 

5U 

\ 

/ 

ZG  ZG 

rn 

UJ 

ll 

r  i  ^ 

K  =P 

r  i  ] 

AGZ(tp)'  =A  GZ((p) 

f  1  1 

U/2  pL\ 

ll/2pisJ 

(l/2pi3F2J 

C.  FORCE  AND  MOMENT  REPRESENTATION 


Using  the  method  described  by  Abkowitz  and  Strom-Tejsen  [Ref.  4],  the  sway 


force,  roll  moment,  and  yaw  moment  can  be  expressed  as  Equations  (2.4). 

Surge  Force:  X  =  f(u,v,r,p,S )  (2.4.a) 

Sway  Force:  Y  =  f2(u,v,r,u,v,r,(p,<p,S^  (2.4. b) 

Roll  Moment:  K  =  f3(u,v,r,  u,v,  r,  </>,  <p,  ip,  c> )  (2.4.c) 

Yaw  Moment:  N  =  f4(u,v,r,u,v,r,<p,<p,S^j  (2.4.d) 

Using  a  third  order  Taylor  expansion  for  f,  f>,  and  f3,  equations  (2.5)  are  obtained. 

f\  =  Xm,u~  +  xvryr  +  4v2  +  XrrK  +  x wfr  +  X  (d)  +  Xo  (2.5. a) 

f2  =  Yvv  +  Yvv  +  Ym,v3  +  Yrr  +  Yrr/  +  Ywrv2r  +  YJ  + 

V  +  yv*vV  +  Y^vt2  +  Yn./<p  +  Y(8)  +  Y0 
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f,  =  Kv  +  X„y  +  K,yr  +  y  r  +  Krry  +  Ayt  +  Kj  +  V  + 

+  K^vf-  +  K^rf  +  K(S)+K0-&GZ  (*) 

ft  =  N„v  +  N,v + N„y  +  N„vr2  +Nrr  +  Ntr  +  Nnrf  +Nf0+  Ay/' + 

N„/f+NxMvf  +JV^V+JV,#y2  +  JV(<?)  +  JV0 


The  coefficients  in  Equation  (2.5)  are  called  the  hydrodynamic  coefficients  and  represent 
the  partial  derivative  of  the  force  or  moment  with  respect  to  the  subscripted  variable.  For 


dY  d  ( d  Y 

example,  Yv  represents  —  and  Yvrr  represents  —  — - 

dv  3vl  dr 


—  — -  .  The  terms  Y0,  K0,  and  N0  are 


the  sway  force,  roll  moment,  and  yaw  moment  induced  by  the  propeller,  respectively. 
Similarly,  Y(S),  N(8),  and  K(S)  represent  the  force  and  moments  induced  by  the  rudder. 
The  righting  moment  due  to  the  static  stability  of  the  ship  is  represented  by  the  term 
A  GZ  {(/)). 

D.  RUDDER  FORCE  AND  MOMENT  REPRESENTATION 


The  expressions  used  to  determine  the  rudder  force  and  moments  are  shown  in 
Equations  (2.6)  and  (2.7)  [Ref.  5].  The  parameters  are  defined  in  Table  3. 


*(£)  =  c*^sin(<S) 

(2. 6. a) 

Y  (A)  =  -(\  +  aH)FN  cos  (A) 

(2.6.b) 

K  ( S)  =  (1  +  aH )  zrFn  cos  ( S) 

(2.6.c) 

N(S)  =  -(xR+aHxH)FNcos(S) 

(2.6.d) 

(2-7. a) 
(2.7  .b) 
(2.7.C) 

(2-7. d) 
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(1- Wp)  +  r|(v  +  xp/ 

vR=rv  +  cSl.r  +  cSrn.r 
f 

aR  =  S  +  taiT1 


')2+Cp*V  +  Cprr} 

(2-7. e) 

3 + w2v 

(2.7  .f) 

\ 

VR 

(2-7. g) 

Table  3  Rudder  Force  and  Moment  Parameters 


Parameter 

Description 

Rudder  to  hull  interaction  coefficient 

1  fn 

Normal  force  action  on  the  rudder 

zr 

z  coordinate  of  point  on  which  rudder  force  Yg  acts 

Xr 

x  coordinate  of  point  on  which  rudder  force  Yg  acts 

XH 

x  coordinate  of  point  on  which  normal  force  FN  acts 

1  A 

Rudder  aspect  ratio 

Ar 

Rudder  area 

1  Vr 

Effective  rudder  inflow  velocity 

aR 

Effective  rudder  inflow  angle 

!  Ur,  Vr 

Components  of  rudder  effective  inflow  velocity 

8 

Constant  in  Equation  (2.7c) 

Up 

Effective  propeller  inflow  velocity 

i  k 

Constant  in  Equation  (2.7c) 

Kx 

Thrust  coefficient 

J 

Advance  coefficient 

n 

Number  of  propeller  revolutions  per  second 

1  D 

Propeller  diameter 

Wp 

Effective  propeller  wake  fraction 

i 

Constant  in  Equation  (2.7e) 

Xp 

x  coordinate  of  propeller  position 

Cpv?  Cpr 

Propeller  flow  rectification  coefficients 

1  Y 

Flow  rectification  coefficient 

Cgrr?  Csrn-,  Cgn-y 

Rudder  wake  coefficients 

E.  FORCES  AND  MOMENTS  FOR  CONTAINER  SHIP  WITH  RUDDER 

The  following  equations  are  used  to  detennine  the  forces  and  moments  acting  on 
the  ship,  including  the  forces  and  moments  caused  by  the  rudder.  Equations  (2.8)  are  used 
in  the  mathematical  model  and  they  combine  all  of  the  above  equations  with  the  rudder 
[Ref.  6],  The  values  of  the  hydrodynamic  derivatives  and  coefficients  are  given  in  Table 
4. 

x=xy+(i -,)T+x,rvr+xy+xy 

Surge  force:  ,  .  (2. 8. a) 

+x^2  +  CrxFn  sin  (<?)  +  (m  +  my )  vr 

Y  =  Yvv  +  Yrr  +  Ypp  +  Y/  +  Yvvvv3  +  Yrr/  +  Ywrv2  r 
Sway  Force:  +Ym.vr 2  +  Yvv/-</>  +  Y^vcjr  +  Yrr/<p  +  Y^rtf  +  (2.8.b) 

(\+aH)  Fn  cos  ( A)  —  ( in  +  mx )  ur 
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(2.8.c) 


K  =  Kvv  +  Krr  +  KpP  +  KJ  +  Kmv 3  +  Kn/  + 

Roll  Moment:  Kvvrv2r  +  Kwrvr2  +  Kvv/-(/>  +  K^vcjr  +  Krr/~(j)  + 

Kr^2  “  ( 1  +  « //  C0S  (  ^)  +  ™JxUr  ~W(GM)0 

N  =  Nvv  +  Nrr  +  Npp  +  N/  +  Nmv 3  +  Nrn /  + 

Yaw  Moment:  Nvvrv2r  +  Nm.vr2  +  N^v2#  +  Nv^v</>2  +  (2.8.d) 

Nrr/2(t>  +  Nr^r</>2  +  (xR  +  aHxH  )Fn  cos(S) 

F.  POD  FORCE  AND  MOMENT  REPRESENTATION 

The  expressions  used  to  detennine  the  pod  force  and  moments  are  shown  in 
Equations  (2.9).  These  equations  were  derived  by  comparing  the  thrust  term  in  the 
equations  of  motion  to  the  rudder  term.  The  thrust  tenns  were  then  modified  by  the  same 
sign,  distance  (for  moments),  and  trigonometric  relationship  as  the  rudder  term,  and  the 
rudder  tenn  was  removed  from  the  equations. 


X„=(l-t)Tsm(Sp^) 

(2-9. a) 

y„=(i-<)rcos(<?^) 

(2-9. b) 

K„=-zr(  l-t)T  COS  (S^) 

(2.9.c) 

(2-9. d) 

G.  FORCES  AND  MOMENTS  FOR  CONTAINER  SHIP  WITH  POD 

The  following  equations  are  used  to  detennine  the  forces  and  moments  acting  on 
the  ship,  including  the  forces  and  moments  caused  by  the  pod.  Equations  (2.10)  are  used 
in  the  mathematical  model  and  they  combine  all  of  the  above  equations  with  the  pod.  The 
values  of  the  hydrodynamic  derivatives  and  coefficients  are  given  in  Table  4. 

x  =  xy  +  (1  ■ - 1 ) T, sin  (S^  )  ■ +  Xwvr + Xwv>  + 

Surge  Force: 

Xm2 +XH>(/)2  +[m  +  my)vr 

Y  =  {\-t)Tcos(8pod)  +  Yvv  +  Yrr  +  YpP  +  YJ  +  Ymv2  + 

Sway  Force:  Yrr/  +  Yvwv2r  +  Yvrrvr 2  +  7^v  V  +  Yv^v<f>2  + 

Yrr/$  +  Y,Hy  ~{m  +  mx )  ur 


(2. 10. a) 


(2. 1 O.b) 
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K  =  -z j 

Roll  Moment: 

N  - 

Yaw  Moment: 


■  (1 "  t)  Tcos  (8pod  )  +  +  Krr  +  KpP  +  KJ  +  Kmv3  + 

Ky  +  Kwrv2r  +  Kwrvr2  +  Kw/</>  +  K^v<t>2  +  (2.10.c) 

Krr/<!)  +  Kr^r<f>2  +  mjxur  -W(GM)(/> 

-S{l-t)TCos(Spod)  +  Nvv  +  Nrr  +  NpP  +  N/  + 

Kvy+NiT/+Nn,rv2r  +  NviTvr2  +  (2.10.d) 

Nw/</>  +  NvHv(/>2  +  Nrr/</>  + 
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Table  4  Hydrodynamic  Derivatives  and  Coefficients  [Ref.  3] 


A.  Hull  only 


i  Yv 

-0.0116000 

Nv 

-0.0038545 

Kv 

0.00030260 

Yr 

0.0024200 

Nr 

-0.0022200 

Kr 

-0.00006300 

Y» 

-0.0000630 

Nj, 

-0.0001424 

K  * 

-0.00002100 

Y 

1  vvv 

-0.1090000 

Nyvv 

0.0014920 

KyvV 

0.00284300 

Y 

i  vrr 

-0.0405000 

NVrr 

0.0015600 

KVrr 

0.00105650 

Y 

1  rrr 

0.0017700 

Nrn- 

-0.0022900 

Km 

-0.00004620 

Y 

i  vvr 

0.0214000 

Nwr 

-0.424000 

Kyyr 

-0.00055800 

YVV(J) 

0.0460500 

Nyv<|) 

-0.0190580 

KVvd) 

-0.00120120 

YV(()(j) 

0.0030400 

Nvdxj) 

-0.0053766 

KV(t)(J) 

-0.00007930 

Yrr(|) 

0.0093250 

Nrrd) 

-0.0038592 

Krrd) 

-0.00024300 

Yr(|)(|) 

-0.0013560 

Nrd)d) 

0.0024195 

0.00003569 

xuu 

-0.0004226 

Xvr 

-0.00311 

Xvv 

-0.00386 

xn 

0.00020 

xw 

-0.00020 

xg 

0.0 

m 

0.00792 

mx 

0.000238 

mv 

0.000238 

i  lx 

0.0000176 

Iv 

0.0 

Iz 

0.000456 

,  Jx 

0.0000034 

CXy 

0.05 

Jz 

0.000419 

lx 

0.0313 

lv 

0.0313 

B.  Propeller  and  rudder 


i  an 

0.237 

T 

1.090 

(1-t) 

0.825 

,  xH 

-0.480 

8 

0.921 

Xr 

-0.500 

!  Crx 

0.710 

k 

0.631 

Xp 

-0.526 

Zr 

0.033 

^5r 

-0.156 

Kj 

0. 527-0. 455J 

CRr 

-0.156 

CRrrr 

-0.275 

CRrrv 

1.960 

CpV 

0.000 

Cgirr 

-0.275 

y 

0.088  v>0 
0.193  v<0 

Cpr 

0.000 

C§rrv 

1.960 

NP 

79.10  Fn=0.2 
118.64  Fn=0.3 
158.19  Fn=0.4 
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H.  DIMENSIONAL  STATE  DERIVATIVES 


Equation  (2.11),  which  defines  the  derivatives  of  the  desired  variables,  is  used  to 
find  the  variation  of  the  following  parameters  over  time:  surge  velocity  u,  sway  velocity 
v,  yaw  velocity  r,  pitch  rate  p,  roll  angle  (j),  heading  angle  vj/,  position  in  x  direction, 
position  in  y  direction,  rudder  deflection  angle  5  (or  5pod),  and  n  (the  actual  shaft 
velocity).  Using  Euler  Integration  to  find  the  desired  variable,  the  derivative  is  multiplied 
by  a  specified  time  step  and  then  added  to  the  previous  (or  initial)  value  of  the  variable. 
The  variable  is  then  equal  to  the  new  calculated  value.  Only  u,  U,  v,  r,  (j),  and  x  and  y 
were  used  for  analysis.  The  equation  is  in  dimensional  form. 


u 

v 

P 

r 

4> 

V 

x 

y 

8 

h 


X 


v_ 

L 


2  A 


/( m  +  mx ) 


( K  +  J,  )(I,  +  J,)K  +  (mray )(/,  +  /,)  N  f  V2\ 


(m  +  mv )(/x  +  Jx )(/_,+  J, ) - (-/«,/, )"  (/,  +J:)~ (m, .a,. )  (lx  +  JX)VL  J 
(  '"J.  )(I2+-J2)Y  +  K  (m  +  my )  (/.  +JZ)-K  (myay  )"  +  (~myly  )(myay )  N  f  v2  ) 


(m  +  my  )(IX  +  Jx  )(/z  +  J, )  -(-/«,/,  )"  (/-  +  J._ )  -(/««,  )  (4  +JX) 

(myay  )  (4  +  Jx  )  Y  +  (~myly  )('«,«,  )  K  +  N  (/«  +  My  )  (I x  +  Jx  )  ~  N(~myly  )7p2) 


(m  +  iiiy  )(IX  +  J x )(/.  +  Jz )  (  /»,./, )“  (/-  +  J-_ )  ~(myay  )  (lx  +  Jx ) 

(V\ 

p  7 


yL2  j 


cos 


\Lj 


(cos(^)w  -sin(^-)cos(^)v)F 
(sin(^-)i/  +  cos(^)cos(^)v)F 


c  _  c 

^ commanded 


'commanded 


-h)x6  0 


(2.11) 


The  commanded  rudder/pod  angle  is  the  desired  deflection  angle  of  the  rudder/pod  for 
the  turn.  The  actual  angle  is  initially  zero  since  the  ship  is  traveling  in  a  straight  line,  and 
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will  increase  each  time  step  as  8  is  integrated.  The  results  obtained  using  Euler 
Integration  are  then  plotted  to  observe  and  compare  the  ship  response  to  a  rudder  versus  a 
podded  propulsor.  The  parameters  compared  are  the  total  speed  over  time,  the  sway 
velocity  over  time,  the  roll  angle  over  time,  and  the  ship’s  position  in  the  xy-plane. 
Analysis  of  the  ship’s  position  in  the  xy-plane  enables  a  comparison  of  the  turning  path 
of  each  ship  configuration.  The  definitions  for  the  turning  path  of  a  ship  are  shown  in 
Figure  2. 


Figure  2.  Turning  Path  of  a  Ship  [Ref.  7] 
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III.  TRANSIENT  ANALYSIS 


A.  INPUT  PARAMETERS 


The  ship’s  behavior  is  simulated  at  various  initial  speeds  U  and  the  rudder  (or 
pod)  angles  as  shown  in  Table  5. 

Table  5  Values  of  Varied  Parameters 


Speed  U  [knots] 

Rudder  or  Pod  Angle  5  [deg] 

5 

5 

10 

10 

15 

15 

20 

20 

1  25 

30 

At  each  speed,  the  rudder  and  pod  angles  are  set  to  the  values  shown  above.  The 
results  are  then  plotted.  The  following  parameters  are  plotted:  total  speed  U  versus  time, 
roll  angle  (j)  versus  time,  sway  velocity  v  versus  time,  and  position  in  the  xy-plane.  Each 
parameter  is  then  studied  at  the  varying  speeds  and  rudder/pod  angles. 

B.  TOTAL  SPEED  U  VERSUS  TIME 

For  the  ship  with  a  pod,  the  total  speed  U  increases  as  pod  angle  increases  (with 
the  initial  speed  constant).  The  speed  U  does  not  vary  as  initial  speed  increases  (with 
constant  pod  angle).  For  the  ship  with  a  rudder,  the  total  speed  U  decreases  as  rudder 
angle  increases  (with  constant  initial  speed).  The  speed  U  again  does  not  change  as  initial 
speed  increases  (with  constant  rudder  angle).  The  behavior  of  the  speed  for  each  ship  is 
shown  in  Table  6  for  each  speed  and  deflection  angle  modeled.  For  both  ships  modeled 
the  speed  of  the  ship  increases  from  its  initial  value  for  an  initial  speed  of  5  knots  and  all 
deflections.  At  an  initial  speed  of  10  and  15  knots,  the  ship  with  the  pod  slows  down  at 
deflections  of  5  and  10  degrees,  while  the  ship  with  the  rudder  speeds  up  (except  at  15 
knots  and  10  degrees  where  it  slows  down).  At  all  other  initial  speeds  and  deflections, 
both  ships  will  slow  down. 
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The  model  indicates  that  for  the  same  initial  speed  U,  the  ship  with  a  rudder  tends 
to  slow  down  as  it  maneuvers  along  a  turn  whereas  the  ship  with  a  pod  tends  to  speed  up. 
This  pod  behavior  is  consistent  with  other  studies  performed  on  the  high  speed 
performance  of  ships  with  pods  installed  [Ref.  1].  The  plots  for  a  total  speed  of  5  knots 
and  varying  rudder/pod  angles  are  shown  in  Figures  3  through  6. 
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Table  6  Ship  Speed  Behavior  at  Different  Initial  Speeds  and 

Deflection  Angles 


Total  Speed  U 
[knots] 

Rudder  or  Pod 
Angle  [deg] 

Ship  Speed  Behavior 
for  Rudder 

Ship  Speed 
Behavior  for  Pod 

5 

5 

t 

t 

5 

10 

t 

t 

5 

15 

t 

T 

5 

20 

t 

t 

10 

5 

T 

1 

10 

10 

T 

1 

10 

15 

t 

t 

10 

20 

t 

t 

15 

5 

t 

1 

15 

10 

1 

1 

!  15 

15 

I 

I 

15 

20 

1 

1 

20 

5 

1 

1 

20 

10 

I 

I 

20 

15 

1 

1 

20 

20 

1 

1 

25 

5 

1 

I 

25 

10 

1 

1 

25 

15 

1 

1 

25 

20 

1 

1 

30 

5 

I 

I 

|  30 

10 

1 

1 

30 

15 

1 

1 

30 

20 

1 

1 
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U  [knots]  U  [knots] 


Total  speed  U 
U  =  5  knots  and  6=5  deg 


Figure  3.  U  vs.  Time,  U  =  5  knots  and  5  =  5° 


Total  speed  U 
U  =  5  knots  and  6=  10  deg 


Figure  4.  U  vs.  Time,  U  =  5  knots  and  8=10° 
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U  [knots]  U  [knots] 


Total  speed  U 
U  =  5  knots  and  6=15  deg 


Figure  5.  U  vs.  Time,  U  =  5  knots  and  5=15° 


Total  speed  U 
U  =  5  knots  and  6  =  20  deg 


Figure  6.  U  vs.  Time,  U  =  5  knots  and  8  =  20° 
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C.  ROLL  ANGLE  4>  VERSUS  TIME 

For  the  ship  with  a  pod,  the  roll  angle  cj)  decreases  as  the  pod  angle  5  increases 
(with  constant  initial  speed  U).  The  roll  angle  cj)  is  approximately  constant  for  varying 
total  speed  U  and  same  pod  angle  5.  For  the  ship  with  a  rudder,  the  results  are  the  same  as 
for  the  pod,  although  the  magnitude  of  the  roll  angle  is  less  than  that  of  the  ship  with  the 
pod  for  the  same  rudder/pod  angle.  The  model  therefore  indicates  that  the  ship  with  the 
pod  will  tend  to  have  a  larger  magnitude  roll  angle  than  the  ship  with  the  rudder.  This 
means  that  the  ship  with  the  pod  will  be  tilted  at  a  steeper  angle  during  a  turn  than  the 
ship  with  the  rudder.  These  results  are  consistent  with  other  studies  performed  on  ships 
with  rudders  [Ref.  1].  The  plots  for  a  total  speed  of  5  knots  and  varying  rudder/pod  angles 
are  shown  in  Figures  7  through  10.  The  plots  for  a  rudder/pod  angle  of  5  degrees  and 
varying  total  speed  are  shown  in  Figures  7  through  16. 
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4>[deg]  <t>[deg] 


Roll  Angle  4> 

U  =  5  knots  and  6=5  deg 


Figure  7.  cj)  vs.  Time,  U  =  5  knots  and  5  =  5° 


Roll  Angle  4> 

U  =  5  knots  and  6=10  deg 


Figure  8.  cj)  vs.  Time,  U  =  5  knots  and  5  =  10° 
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<t>  [deg]  4,  [deg] 


4 


Roll  Angle  4> 

U  =  5  knots  and  6=15  deg 


Figure  9.  cj)  vs.  Time,  U  =  5  knots  and  5=15° 


Roll  Angle  4> 

U  =  5  knots  and  6  =  20  deg 


Figure  10.  cj)  vs.  Time,  U  =  5  knots  and  5  =  20° 
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Roll  Angle  <f> 

U  =  5  knots  and  6  =  5  deg 


Figure  11.  cj)  vs.  Time,  U  =  5  knots  and  8  =  5°  (same  as  Figure  7) 


Roll  Angle  4> 

U  =  10  knots  and  8=5  deg 


Figure  12.  (j)  vs.  Time,  U  =  10  knots  and  5  =  5° 
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<(>  [deg]  4,  [deg] 


Roll  Angle  4> 

U  =  15  knots  and  5=5  deg 


Figure  13.  cj)  vs.  Time,  U  =  15  knots  and  5  =  5° 


Roll  Angle  4> 

U  =  20  knots  and  6  =  5deg 


Figure  14.  cj)  vs.  Time,  U  =  20  knots  and  5  =  5° 
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4 


Roll  Angle  4> 

U  =  25  knots  and  8=5  deg 


Figure  15.  cj)  vs.  Time,  U  =  25  knots  and  5  =  5° 


Roll  Angle  4> 

U  =  30  knots  and  8=5  deg 


Figure  16.  cj)  vs.  Time,  U  =  30  knots  and  5  =  5° 
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D.  SWAY  VELOCITY  V  VERSUS  TIME 

For  the  ship  with  a  pod,  the  sway  velocity  v  is  approximately  constant  as  the 
speed  U  increases  (with  constant  pod  angle  5).  The  sway  velocity  v  increases  in 
magnitude  as  the  pod  angle  5  increases  (with  constant  speed  U).  For  the  ship  with  a 
rudder,  the  sway  velocity  v  is  constant  as  the  initial  speed  U  increases  (with  constant 
rudder  angle  5).  The  sway  velocity  v  increases  in  magnitude  as  the  rudder  angle  8 
increases  (with  constant  speed  U).  The  plots  for  a  speed  of  5  knots  and  varying 
rudder/pod  angle  are  shown  in  Figures  17  through  20. 

As  seen  from  the  plots,  the  sway  velocity  is  always  negative  for  the  rudder  and  is 
negative  for  pod  angles  greater  than  5°,  regardless  of  the  speed.  However,  the  magnitude 
of  the  sway  velocity  is  always  larger  for  the  rudder  than  for  the  pod. 
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v  [knots]  v  [knots] 


Sway  Velocity  v 


Figure  17.  v  vs.  time,  U  =  5  knots  and  5  =  5° 


Sway  Velocity  v 
U  =  5  knots  and  5=  10  deg 


Figure  18.  v  vs.  time,  U  =  5  knots  and  5  =  10° 
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v  [knots]  v  [knots] 


1 


Sway  Velocity  v 
U  =  5  knots  and  6=  15  deg 


Figure  19.  v  vs.  time,  U  =  5  knots  and  5=15° 


Sway  Velocity  v 
U  =  5  knots  and  6  =  20  deg 


Figure  20.  v  vs.  time,  U  =  5  knots  and  5  =  20° 
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E. 


POSITION  IN  XY  PLANE  VERSUS  TIME 


The  ship  with  a  pod  has  a  smaller  turning  radius  than  the  ship  with  a  rudder.  The 
ship  with  the  pod  also  has  a  smaller  advance  and  transfer  at  the  commanded  pod  angle  for 
increasing  pod  angles  until  the  total  speed  is  15  knots  and  the  rudder/pod  angle  is  15 
degrees.  After  this  point,  the  ship  with  the  rudder  has  a  smaller  advance.  Also,  as  the 
rudder/pod  angle  increases,  the  pod  progressively  has  a  larger  turning  radius  than  before 
for  the  same  rudder/pod  angle  (and  different  speed  U).  However,  the  ship  with  the  pod 
still  has  a  smaller  turning  radius  than  the  ship  with  a  rudder.  As  the  total  speed  U 
increases,  the  advance  of  the  ship  with  the  pod  is  greatly  increased.  The  advance 
increases  to  as  much  as  2000  m  for  a  total  speed  of  30  knots  and  pod  angle  of  20  degrees. 
The  transfer  also  increases  so  that  at  higher  speeds  U  the  pod  requires  more  distance  to 
begin  the  turn.  These  results  are  consistent  with  previous  studies  [Ref.  1].  The  plots  for 
varying  rudder/pod  angle  and  varying  total  speed  U  are  shown  in  Figures  21  through  32. 
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Position  in  the  xy  Plane 
U  =  5  knots  and  6=5  deg 


Figure  2 1 .  Position  in  XY  Plane,  U  =  5  knots  and  5  =  5° 


Position  in  the  xy  Plane 
U  =  5  knots  and  6  =  10  deg 


Figure  22.  Position  in  XY  Plane,  U  =  5  knots  and  5=10° 
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Y  [m]  Y  [m] 


Position  in  the  xy  Plane 
U  =  5  knots  and  6=  15  deg 


Figure  23.  Position  in  XY  Plane,  U  =  5  knots  and  5=15° 


Position  in  the  xy  Plane 
U  =  5  knots  and  5  =  20  deg 


Figure  24.  Position  in  XY  Plane,  U  =  5  knots  and  5  =  20° 
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Y  [m]  Y  [m] 


2500 


Position  in  the  xy  Plane 
U  =  15  knots  and  5=5  deg 


Figure  25.  Position  in  XY  Plane,  U  =  15  knots  and  5  =  5° 


Position  in  the  xy  Plane 
U  =  15  knots  and  6=  10  deg 


Figure  26.  Position  in  XY  Plane,  U  =  15  knots  and  5=10° 
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Y  [m]  Y  [m] 


Position  in  the  xy  Plane 
U  =  15  knots  and  6=15  deg 


Figure  27.  Position  in  XY  Plane,  U  =  15  knots  and  5=15° 


Position  in  the  xy  Plane 
U  =  15  knots  and  6  =  20  deg 


Figure  28.  Position  in  XY  Plane,  U  =  15  knots  and  5  =  20° 
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Y  [m]  Y  [m] 


Position  in  the  xy  Plane 
U  =  30  knots  and  5  =  5  deg 


Figure  29.  Position  in  XY  Plane,  U  =  30  knots  and  5  =  5° 


Position  in  the  xy  Plane 
U  =  30  knots  and  6=  10  deg 


Figure  30.  Position  in  XY  Plane,  U  =  30  knots  and  5=10° 
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Y  [m]  Y  [m] 


1400 


Position  in  the  xy  Plane 
U  =  30  knots  and  6=15  deg 


Figure  3 1 .  Position  in  XY  Plane,  U  =  30  knots  and  5=15° 


Position  in  the  xy  Plane 
U  =  30  knots  and  6  =  20  deg 


Figure  32.  Position  in  XY  Plane,  U  =  30  knots  and  5  =  20° 
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IV.  STEADY  STATE  ANALYSIS 


A.  PARAMETERS  OF  INTEREST 

The  final  integrated  value  of  each  parameter  found  from  numerical  integration  of 
Equation  2.11  is  obtained  from  Matlab;  this  value  is  the  steady  state  value  of  the 
parameter.  For  the  steady  state  analysis,  the  following  parameters  are  studied:  sway 
velocity  v,  yaw  velocity  r,  and  roll  angle  cj). 

B.  ANALYSIS  OF  PARAMETERS 

The  results  are  plotted  for  different  initial  speeds  U  and  different  rudder/pod 
angles  5.  As  observed  from  the  plots,  the  steady  state  values  of  each  parameter  do  not 
vary  with  a  change  in  the  initial  speed  U  (with  the  same  rudder/pod  angle).  This  is  the 
same  result  obtained  from  the  transient  analysis  for  the  parameters  studied.  The 
parameters  vary  slightly  when  plotted  using  constant  initial  speed  U  and  varying 
rudder/pod  angle. 
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Steady  State  Sway  Velocity  v 


Figure  33. 


Steady  State  Sway  Velocity  for  Pod  with  Varying  Initial  Speed 
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Pod  Angle  [deg] 


Steady  State  Sway  Velocity  v 


Figure  34. 


Steady  State  Sway  Velocity  for  Rudder  with  Varying  Initial  Speed 
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Rudder  Angle  [deg] 


Steady  State  YawVelocityr 


Figure  35. 


Steady  State  Yaw  Velocity  for  Pod  with  Varying  Initial  Speed 
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Pod  Angle  [deg] 


Steady  State  Yaw  Velocity  r 


Figure  36. 


Steady  State  Yaw  Velocity  for  Rudder  with  Varying  Initial  Speed 
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Rudder  Angle  [deg] 


Steady  State  Roll  Angle  phi 


Figure  37. 


Steady  State  Roll  Angle  for  Pod  with  Varying  Initial  Speed 
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Pod  Angle  [deg] 


Steady  State  Roll  Angle  phi 


Figure  38. 


Steady  State  Roll  Angle  for  Rudder  with  Varying  Initial  Speed 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

In  this  thesis,  a  mathematical  model  was  modified  to  simulate  the  motion  of  a  ship 
operating  with  a  rudder  and  operating  with  a  pod.  The  model  was  first  used  to  simulate 
the  ship  operating  with  a  rudder.  The  model  was  then  modified  to  simulate  operation  with 
a  pod.  The  two  results  were  then  compared  to  detennine  how  well  a  ship  with  a  pod 
operates  at  various  speeds  and  pod  deflections.  Based  on  the  comparison  made,  the 
following  results  were  observed  for  the  ship  with  the  pod: 

1.  The  results  suggest  that  a  ship  with  a  pod  will  operate  very  well  at  speeds  less 
than  or  equal  to  15  knots  and  pod  deflections  of  less  than  or  equal  to  15  degrees. 
After  this  point,  the  pod  tends  to  have  a  larger  advance  and  transfer  as  initial 
speed  and  deflection  angle  increase. 

2.  The  results  indicate  that  for  a  ship  with  a  pod,  the  ship’s  total  speed  will  increase 
during  a  turn  if  the  initial  speed  is  5  knots  or  less.  At  all  other  speeds  and 
deflections,  the  ship  with  the  pod  will  slow  down.  This  behavior  is  similar  to  that 
of  a  ship  with  a  rudder. 

3.  The  results  show  that  a  ship  with  a  pod  will  tend  to  have  a  larger  magnitude  roll 
angle  than  a  ship  with  a  rudder;  however,  the  roll  angle  does  not  exceed  8° 
regardless  of  the  initial  speed. 

4.  The  results  indicate  that  the  sway  velocity  v  will  be  approximately  constant  as  the 
initial  speed  U  increases  if  the  pod  deflection  is  constant.  The  sway  velocity  will 
decrease  as  the  pod  deflection  increases  and  the  initial  speed  is  held  constant. 

Based  on  the  results  observed  using  the  model  as  described,  the  rudder  performs  better 
than  the  pod  at  high  speeds  and  high  deflection  angles.  Naval  vessels  must  have  the 
capability  of  responding  rapidly  in  an  emergency,  regardless  of  the  initial  speed  of  the 
vessel.  The  pod  does  not  provide  the  necessary  performance  requirements  to  meet  this 
constraint  in  the  current  configuration. 
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B.  RECOMMENDATIONS 

The  results  produced  in  this  thesis  were  intended  as  an  initial  study  into  the 
viability  of  replacing  a  rudder  on  a  naval  ship  with  a  podded  propulsor.  In  order  to 
accomplish  this,  an  existing  model  was  modified  to  simulate  a  podded  propulsor  in  place 
of  a  propeller  and  no  rudder.  The  model  did  not  take  into  account  any  change  in  the  hull 
configuration  that  may  be  required  for  use  of  a  pod,  nor  the  possible  need  of  more  than 
one  pod  to  properly  maneuver  the  ship.  One  area  of  further  study  is  to  further  modify  the 
model  used  to  account  for  this,  or  to  create  a  new  model  to  simulate  a  ship  that  already 
uses  a  pod. 

Further  study  is  also  necessary  to  determine  the  best  possible  pod  and  hull 
configuration  for  a  given  application.  Previous  studies  have  shown  that  pods  do  not 
operate  as  well  as  rudders  at  high  speeds.  Further  study  is  needed  to  investigate  the 
benefits  of  adding  control  surfaces  to  counteract  this  negative  perfonnance  aspect  of  a 
pod. 


48 


APPENDIX  A.  MATLAB  CODE  FOR  SHIP  WITH  RUDDER 


A.  V07CONT 

function  ul  =  V07_cont(k,x,u) 
ul  =  [  20*pi/180;  %  delta  (rad) 

85  ]; 

B.  V07DATA 

%  V07  DATA  datafile  for  v07.m.  Contains  initial  state 
%  vector,  initial  input  and  main  dimensions. 

% 

%  NTH  1994  Trygve  Lauvdal 
clc 

echo  on 

%  Default  parameters  for  High  Speed  Container  Ship: 

% 

%  Max  rudder  angle  :  0  deg 
%  Max  rudder  rate  :  0  deg/s 
%  Max  shaft  velocity  :  160  rpm 
%  Length  of  ship  :  175  m 

%  Initial  states,  inputs  and  sample  time  (to  be  edited  by  user): 


u  = 

15.433; 

%  surge  velocity 

(m/s) 

V  = 

0; 

%  sway  velocity 

(m/s) 

p 

0; 

%  roll  velocity 

(rad/s) 

r  = 

0; 

%  yaw  velocity 

(rad/s) 

phi  = 

0; 

%  roll  angle 

(rad) 

psi  = 

0; 

%  yaw  angle 

(rad) 

xpos  = 

0; 

%  surge  position 

(m) 

ypos  = 

0; 

%  sway  position 

(m) 

delta  = 

0; 

%  actual  rudder  angle  (rad) 

n  = 

10; 

%  actual  shaft  speed 

(rpm) 
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d_c 

n_c 

h 


20*pi/180; 

0; 

0.1; 


%  commanded  rudder  angle  (rad) 
%  commanded  shaft  speed  (rpm) 

%  sample  time  (s) 


x 

u 


[u  v  p  r  phi  psi  xpos  ypos  delta  n]'; 
[d_c  n_c]'; 


%  state  vector 
%  input  vector 


echo 

current  =  pwd; 

if  ~exist(['sim.m’]), 

saveon  =  [current  7v07_data  xuh']; 

else 

saveon  =  [current  7v07/v07_data  xuh’]; 
end 

eval(['save  v07_data’]) 
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C.  V07RUDDER 

function  xdot  =  V07_rudder(x,ui) 

% 

%  xdot  =  V07(x,ui);  returns  the  time  derivate  of  the  state  vector 

% 

%  x  =  [  u  v  p  r  phi  psi  xpos  ypos  delta  n  ]’  where 

% 


%u  = 

surge  velocity 

(m/s) 

%v  = 

sway  velocity 

(m/s) 

%p  = 

roll  velocity 

(rad/s) 

%  r  = 

yaw  velocity 

(rad/s) 

%  phi  = 

roll  angle 

(rad) 

%  psi  = 

yaw  angle 

(rad) 

%  xpos  = 

position  in  x-direction 

(m) 

%  ypos  = 

position  in  y-direction 

(m) 

%  delta  = 

actual  rudder  angle 

(rad) 

%  n  = 

actual  shaft  velocity 

(rpm) 

% 


%  The  input  vector  is: 

% 

%  ui  =  [  delta_c  n_c  ]’  where 

% 

%  delta  c  =  commanded  rudder  angle  (rad) 

%  n_c  =  commanded  shaft  velocity  (rpm) 

% 

%  NTH  1994,  Trygve  Lauvdal 

% 

%  Reference  :  Son  og  Nomoto  (1982). 

%  On  the  Coupled  Motion  of  Steering  and 

%  Rolling  of  a  High  Speed  Container  Ship, 

%  Naval  Architect  of  Ocean  Engineering, 

%  20:  73-83.  From  J.S.N.A. ,  Japan,  Vol.  150,  1981. 
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%  Check  of  input  and  state  dimensions 

if  ~(length(x)  ==  10),error('x-vector  must  have  dimension  10  !');end 
if -(length(ui)  ==  2), error('u-vector  must  have  dimension  2  !’);end 

%  Normalization  variables 

L  =  175;  %  length  of  ship  (m) 

sqrV  =  (x(l)A2  +  x(2)A2); 

V  =  sqrt(sqrV);  %  service  speed  (m/s) 

%  Check  of  service  speed 

if  V  ==  0,error(’The  ship  must  have  speed  greater  than  zero’);end 


delta  max  = 

20; 

%  max  rudder  angle 

(deg) 

D  delta  max  = 

5; 

%  max  rudder  angle  rate 

(deg/s) 

n  max  = 

160; 

%  max  shaft  velocity 

(rpm) 

%  Non-dimensional  states  and  inputs 

deltac  =  ui(l); 

n_c  =  ui(2)/60*L/V; 


u  =x(l)/V;  v  =  x(2)/V; 

p  =  x(3)*L/V;r  =  x(4)*L/V; 

phi  =  x(5);  psi  =  x(6); 
delta  =  x(9);  n  =  x(10)/60*L/V; 


%  Parameters,  hydrodynamic  derivatives  and  main  dimensions 
m  =  0.00792;  mx  =  0.000238;  my  =  0.000238; 

lx  =  0.0000176;  alphay  =0.05;  lx  =  0.03 13; 

ly  =  0.0313;  lx  =0.0000176;  Iz  =  0.000456; 

Jx  =  0.0000034;  Jz  =0.000419;  xG  =  0; 


B  =25.40; 


dF  =  8.00;  g 
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=  9.81; 


dA  =9.00;  d  =  8.50;  nabla  =  2 1222;%Displacement  Volume 

KM  =10.39;  KB  =  4.6154;  AR  =  33.0376;%Rudder  Area 
KG  =10.09;  Delta  =  1.8219;  %Aspect  Ratio 

D  =  6.533;  %Propeller  Diameter 
GM  =  ,3/L; 

rho  =  1000;  t  =0.175;  T  =0.0005; 


W  =  rho*g*nabla/(rho*LA2*VA2/2); 


Xuu  =-0.0004226;  Xvr  =-0.00311;  Xrr  =0.00020; 
Xphiphi  =  -0.00020;  Xvv  =  -0.00386; 


Kv  =  0.0003026; 
Kphi  =-0.000021; 
Kvvr  =  -0.000588; 
Kvphiphi  =  -0.0000793; 


Kr  =  -0.000063; 
Kvvv  =  0.002843; 
Kvrr  =  0.0010565; 
Krrphi  = -0.000243; 


Kp  =  -0.0000075; 
Krrr  =  -0.0000462; 
Kvvphi  =-0.0012012; 
Krphiphi  =  0.00003569; 


Yv  =-0.0116;  Yr  =  0.00242; 
Yphi  =-0.000063;  Yvvv  =-0.109; 
Yvvr  =  0.0214;  Yvrr  =-0.0405; 
Yvphiphi=  0.00304;  Yrrphi=  0.009325; 


Yp  =  0; 

Yrrr  =  0.00177; 

Yvvphi  =  0.04605; 
Yrphiphi  =  -0.001368; 


Nv  = -0.0038545;  Nr  =-0.00222;  Np  =  0.000213; 

Nphi  = -0.000 1424;Nvvv  =  0.001492;  Nrrr  =-0.00229; 

Nvvr  =-0.0424;  Nvrr  =  0.00156;  Nvvphi  =-0.019058; 
Nvphiphi  = -0.0053766;  Nrrphi  = -0.0038592;  Nrphiphi  =  0.0024195; 


kk  =  0.631; 
wp  =  0.184; 
cpv  =  0.0; 
cRr  =-0.156; 
cRX  =  0.71; 
xH  =  -0.48; 


epsilon  =  0.921; 
tau  =  1.09; 
cpr  =  0.0; 
cRrrr  =-0.275; 
aH  =  0.237; 


xR  =  -0.5; 
xp  =  -0.526; 
ga  =  0.088; 
cRrrv=  1.96; 
zR  =  0.033; 
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%  Masses  and  moments  of  inertia 

ml  1  =  (m+mx); 

m2  2  =  (m+my); 

m32  =  -my*ly; 

m42  =  my*alphay; 

m3  3  =  (Ix+Jx); 

m44  =  (Iz+Jz); 

%  Rudder  saturation  and  dynamics 
if  abs(delta  c)  >=  delta_max*pi/180, 
delta_c  =  sign(delta_c)*delta_max*pi/180; 
end 

deltadot  =  delta  c  -  delta; 
if  abs(delta  dot)  >=  Ddelta_max*pi/180, 
deltadot  =  s  ign(deltadot)  *  D  delta_max*pi  / 180; 
end 

%  Shaft  velocity  saturation  and  dynamics 
n_c  =  n_c*V/L; 
n  =  n*V/L; 
if  abs(n_c)  >=  n_max/60, 
n_c  =  sign(n_c)*n_max/60; 
end 

ifn  >  0.3,Tm=5.65/n;else,Tm=18.83;end 
ndot  =  l/Tm*(n_c-n)*60; 

%  Calculation  of  state  derivatives 
vR  =  ga*v  +  cRr*r  +  cRrrr*rA3  +  cRrrv*rA2*v; 
uP  =  cos(v)*((l  -  wp)  +  tau*((v  +  xp*r)A2  +  cpv*v  +  cpr*r)); 
J  =  uP*V/(n*D); 

KT  =0.527  -0.455*1; 
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uR  =  uP*epsilon*sqrt(l  +  8*kk*KT/(pi*JA2)); 
alphaR  =  delta  +  atan(vR/uR); 

FN  =  -  ((6.13*Delta)/(Delta  +  2 . 2  5  )  )  *  ( AR/L  A  2  )  *  (uR A 2  +  vRA2)*sin( alphaR); 
T  =  2*rho*DA4/(VA2*LA2*rho)*KT*n*abs(n); 

%  Forces  and  moments 

X  =  Xuu*uA2  +  (l-t)*T  +  Xvr*v*r  +  Xvv*vA2  +  Xrr*rA2  +  Xphiphi*phiA2  +  ... 
(m  +  my)*v*r  +  cRX*FN*sin(  delta); 

Y  =  Yv*v  +  Yr*r  +  Yp*p  +  Yphi*phi  +  Yvvv*vA3  +  Yrrr*rA3  +  ... 

Yvvr*vA2*r  +  Yvrr*v*rA2  +  Yvvphi*vA2*phi  +  Yvphiphi*v*phiA2  + 
Yrrphi*rA2*phi  +  ... 

Yrphiphi*r*phiA2  -  (m  +  mx)*u*r  +  (1  +  aFI)*FN*cos(delta); 

K  =  Kv*v  +  Kr*r  +  Kp*p  +  Kphi*phi  +  Kvvv*vA3  +  Krrr*rA3  +  ... 

Kvvr*vA2*r  +  Kvrr*v*rA2  +  Kvvphi*vA2*phi  +  Kvphiphi*v*phiA2  + 
Krrphi*rA2*phi  +  ... 

Krphiphi*r*phiA2  +  mx*lx*u*r  -  W*GM*phi  -  (1  +  FI)*zR*FN*cos(delta); 

N  =  Nv*v  +  Nr*r  +  Np*p  +  Nphi*phi  +  Nvvv*vA3  +  Nrrr*rA3  +  ... 

Nvvr*vA2*r  +  Nvrr*v*rA2  +  Nvvphi*vA2*phi  +  Nvphiphi*v*phiA2  + 
Nrrphi*rA2*phi  +  ... 

Nrphiphi*r*phiA2  +  (xR  +  aH*xFI)*FN*cos(delta); 

%  Dimensional  state  derivatives 

xdot  =[  X*(VA2/L)/ml  1 

-(-m33*m44*Y+m32*m44*K+m42*m33*N)/(m22*m33*m44-m32A2*m44- 

m42A2*m33)*(VA2/L) 

(-m32*m44*Y+K*m22*m44-K*m42A2+m32*m42*N)/(m22*m33*m44- 
m3  2  A2  *  m44  -m42  A2  *  m3  3  )  *  ( VA2/L  A2) 

(-m42*m33*Y+m32*m42*K+N*m22*m33-N*m32A2)/(m22*m33*m44- 
m3  2  A2  *  m44  -m42  A2  *  m3  3  )  *  ( VA2/L  A2) 

p*(V/L) 

cos(phi)*r*(V/L) 
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(cos(psi)  *u-sin(psi)  *  cos(phi)  *  v)*  V 
(sin(psi)*u+cos(psi)*cos(phi)*v)*V 
deltadot 
ndot  ]; 

D.  V07SIM 

%  SIM.m 

% 

%  Main  program  for  contship  simulation 

% 

%  NTH  1994,  Trygve  Lauvdal 
clear 


v07_data 
loadv07  data 


%  set  model  data 
%  load  model  data 


N  =  15000;  %  number  of  samples 

table  =  zeros(N,length(x)+length(u));%  storage  vector 


disp('Simulating'); 

for  i  =  1:N+1, 
table(i,:)  =  [  x' ,  u'] ; 
u  =  v07_cont(i*h,x,u); 
xdot  =  v07_rudder(x,u'); 
x  =  x  +  h*xdot ; 
end 


%  START  OF  MAIN  LOOP 
%  data  storage 
%  input 

%  xdot(i)  =  f(x(i),u(i)) 

%  euler  integration  x(i+l) 

%  END  OF  MAIN  LOOP 
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E. 


V07RUDDER  PLOT 


function  plotdata  =  V07rudder_plot(table,h,N) 

figure  (1) 
x  =  table(:,7); 
y  =  table(:,8); 
plot(x,y) 

legend('Pod','Rudder') 

figure(2) 
time  =  (0:N)*h; 
u  =  table(:,l)*l. 94384449244; 
v  =  table(:,2)*  1 .94384449244; 

U  =  sqrt(u.*u  +  v.*v); 

plot(time,U), 

legend('Pod','Rudder') 

figure(3) 

plot(time,v), 

legend('Pod','Rudder') 

figure(4) 

phi  =  table(:,5)*(  180/pi); 

plot(time,phi) 

legend('Pod','Rudder') 
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APPENDIX  B.  MATLAB  CODE  FOR  SHIP  WITH  POD 


A.  V07CONT 

function  ul  =  V07_cont(k,x,u) 

ul  =  [  20*pi/180;  %  delta_pod  (rad) 

85  ]; 

B.  V07DATA 

%  V07  DATA  datafile  for  v07.m.  Contains  initial  state 
%  vector,  initial  input  and  main  dimensions. 

% 

%  NTH  1994  Trygve  Lauvdal 
clc 

echo  on 


%  Default  parameters  for  High  Speed  Container  Ship  : 

% 

%  Max  pod  angle 

0  deg 

%  Max  pod  rate 

0  deg/s 

%  Max  shaft  velocity 

160  rpm 

%  Length  of  ship 

175  m 

%  Initial  states,  inputs  and  sample  time  (to  be  edited  by  user): 

u  = 

15.433; 

%  surge  velocity 

(m/s) 

v  = 

0; 

%  sway  velocity 

(m/s) 

P 

0; 

%  roll  velocity 

(rad/s) 

r  = 

0; 

%  yaw  velocity 

(rad/s) 

phi  = 

0; 

%  roll  angle 

(rad) 

psi  = 

0; 

%  yaw  angle 

(rad) 

xpos  = 

0; 

%  surge  position 

(m) 

ypos  = 

0; 

%  sway  position 

(m) 

delta_pod= 

0; 

%  actual  rudder  angle  (rad) 

n  = 

10; 

%  actual  shaft  speed 

(rpm) 

d_c  = 

20*pi/l  80; 

%  commanded  rudder  angle  (rad) 
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n  c 


0; 


%  commanded  shaft  speed  (rpm) 


h  =  0.1;  %  sample  time  (s) 

x  =  [u  v  p  r  phi  psi  xpos  ypos  delta  n]'; 

u  =  [d_c  n_c]’; 

echo 

current  =  pwd; 

if  ~exist(['sim.m’]), 

saveon  =  [current  '/v07_data  xuh']; 

else 

saveon  =  [current  7v07/v07_data xuh’]; 
end 

eval(  ['save  v07_data’]) 


%  state  vector 
%  input  vector 
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c. 


V07  POD 


function  xdot  =  V07_pod(x,ui) 

% 

%  xdot  =  V07(x,ui);  returns  the  time  derivate  of  the  state  vector 

% 


%  x  =  [  u  v  p  r  phi  psi  xpos  ypos  delta  n  ]'  where 

% 


% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 


u  = 

surge  velocity 

(m/s) 

V  = 

sway  velocity 

(m/s) 

p  = 

roll  velocity 

(rad/s) 

r  = 

yaw  velocity 

(rad/s) 

phi  = 

roll  angle 

(rad) 

psi  = 

yaw  angle 

(rad) 

xpos  = 

position  in  x-direction 

(m) 

ypos  = 

position  in  y-direction 

(m) 

delta_pod= 

actual  pod  angle 

(rad) 

n  = 

actual  shaft  velocity 

(rpm) 

The  input  vector  is: 

ui  =  [  deltac  n_c  ]'  where 

deltac  =  commanded  pod  angle 
n_c  =  commanded  shaft  velocity 

NTH  1994,  Trygve  Lauvdal 


(rad) 

(rpm) 


%  Reference  :  Son  og  Nomoto  (1982). 

%  On  the  Coupled  Motion  of  Steering  and 

%  Rolling  of  a  High  Speed  Container  Ship, 

%  Naval  Architect  of  Ocean  Engineering, 

%  20:  73-83.  From  J.S.N.A. ,  Japan,  Vol.  150,  1981. 
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%  Check  of  input  and  state  dimensions 

if  ~(length(x)  ==  10),error('x-vector  must  have  dimension  10  !');end 
if -(length(ui)  ==  2), error('u-vector  must  have  dimension  2  !’);end 

delta_pod  =  20;  %  pod  angle  (deg) 

delta_pod  =  delta_pod*pi/180;  %  pod  angle  (rad) 

%  Normalization  variables 

L  =  175;  %  length  of  ship  (m) 

sqrV  =  (x(l)A2  +  x(2)A2); 

V  =  sqrt(sqrV);  %  service  speed  (m/s) 

%  Check  of  service  speed 

if  V  ==  0,error(’The  ship  must  have  speed  greater  than  zero’);end 


delta  max  = 

20; 

%  max  rudder  angle 

(deg) 

D  delta  max  = 

5; 

%  max  rudder  angle  rate 

(deg/s) 

n  max  = 

160; 

%  max  shaft  velocity 

(rpm) 

%  Non-dimensional  states  and  inputs 

deltac  =  ui(l); 

n_c  =  ui(2)/60*L/V; 

u  =x(l)/V;  v  =x(2)/V; 

p  =  x(3)*L/V;  r  =  x(4)*L/V; 

phi  =  x(5);  psi  =  x(6); 

delta_pod  =  x(9);  n  =  x(10)/60*L/V; 

%  Parameters,  hydrodynamic  derivatives  and  main  dimensions 

m  =  0.00792;  mx  =  0.000238;  my  =  0.000238; 

lx  =  0.0000176;  alphay  =0.05;  lx  =  0.03 13; 

ly  =  0.0313;  lx  =0.0000176;  Iz  =  0.000456; 
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Jx  =  0.0000034;  Jz 


=  0.000419; 


xG  =  0; 


B 

dA 

KM 

KG 

D 

GM 

rho 


=  25.40;  dF  =  8.00;  g  =9.81; 


=  9.00;  d  =  8.50;  nabla  =  2 1222;%Displacement  Volume 

=  10.39;  KB  =  4.6154;  AR  =  33.0376;%Rudder  Area 

=  10.09;  Delta  =  1.8219;  %Aspect  Ratio 
=  6.533;  %Propeller  Diameter 
=  ,3/L; 

=  1000;  t  =0.175;  T  =0.0005; 


W  =  rli  o  *  g  *  n  ab  I  a/(rh  o  *  LA2  *  V  A2/2 ) ; 


Xuu  =-0.0004226;  Xvr  =-0.00311;  Xrr  =0.00020; 
Xphiphi  =  -0.00020;  Xvv  =  -0.00386; 


Kv  =  0.0003026; 
Kphi  =-0.000021; 
Kvvr  =  -0.000588; 
Kvphiphi  =  -0.0000793; 


Kr  =  -0.000063; 
Kvvv  =  0.002843; 
Kvrr  =  0.0010565; 
Krrphi  = -0.000243; 


Kp  =  -0.0000075; 
Krrr  =  -0.0000462; 
Kvvphi  =-0.0012012; 
Krphiphi  =  0.00003569; 


Yv  =-0.0116;  Yr  =  0.00242; 
Yphi  =-0.000063;  Yvvv  =-0.109; 
Yvvr  =  0.0214;  Yvrr  =-0.0405; 
Yvphiphi=  0.00304;  Yrrphi=  0.009325; 


Yp  =  0; 

Yrrr  =  0.00177; 

Yvvphi  =  0.04605; 
Yrphiphi  =  -0.001368; 


Nv  = -0.0038545;  Nr  =-0.00222;  Np  =  0.000213; 

Nphi  = -0.000 1424;Nvvv  =  0.001492;  Nrrr  =-0.00229; 

Nvvr  =-0.0424;  Nvrr  =  0.00156;  Nvvphi  =-0.019058; 
Nvphiphi  = -0.0053766;  Nrrphi  = -0.0038592;  Nrphiphi  =  0.0024195; 


kk 

=  0.631; 

epsilon  =  0.921; 

xR 

=  -0.5; 

wp 

=  0.184; 

tau  =  1.09; 

xp 

=  -0.526; 

cpv 

=  0.0; 

cpr  =  0.0; 

ga 

=  0.088; 
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cRr  =-0.156;  cRrrr  =-0.275;  cRrrv  = 

cRX  =  0.71;  aH  =  0.237;  zR  = 

xH  =  -0.48; 

%  Masses  and  moments  of  inertia 

mil  =  (m+mx); 

m2  2  =  (m+my); 

m3  2  =  -my*ly; 

m42  =  my*alphay; 

m3  3  =  (Ix+Jx); 

m44  =  (Iz+Jz); 

%  Pod  saturation  and  dynamics 
if  abs(delta  c)  >=  delta_max*pi/180, 
delta_c  =  sign(delta_c)*delta_max*pi/180; 
end 

deltadot  =  delta  c  -  delta_pod; 
if  abs(delta  dot)  >=  Ddelta_max*pi/180, 
delta  dot  =  s  ign(deltadot)  *D  delta_max*pi  / 180; 
end 

%  Shaft  velocity  saturation  and  dynamics 
n_c  =  n_c*V/L; 
n  =  n*V/L; 
if  abs(n_c)  >=  n_max/60, 
n_c  =  sign(n_c)*n_max/60; 
end 

if  n  >  0.3,Tm=5.65/n;else,Tm=18.83;end 
ndot  =  l/Tm*(n_c-n)*60; 

%  Calculation  of  state  derivatives 

vR  =  ga*v  +  cRr*r  +  cRrrr*rA3  +  cRrrv*rA2*v; 
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1.96; 

0.033; 


uP  =  cos(v)*((l  -  wp)  +  tau*((v  +  xp*r)A2  +  cpv*v  +  cpr*r)); 

J  =  uP*V/(n*D); 

KT  =  0.527  -0.455*J; 

uR  =  uP*epsilon*sqrt(l  +  8*kk*KT/(pi*JA2)); 

T  =  2*rho*DA4/(VA2*LA2*rho)*KT*n*abs(n); 


%  Forces  and  moments 

X  =  Xuu*uA2  +  ( 1  -t)  *  T  *  sin(delta_pod)  +  Xvr*v*r  +  Xvv*vA2  +  Xrr*rA2  + 
Xphiphi*phiA2  +  (m  +  my)*v*r; 


Y  =  (l-t)*T*cos(delta_pod)  +  Yv*v  +  Yr*r  +  Yp*p  +  Yphi*phi  +  Yvvv*vA3  + 
Yrrr*rA3  +  ... 

Yvvr*vA2*r  +  Yvrr*v*rA2  +  Yvvphi*vA2*phi  +  Yvphiphi*v*phiA2  + 
Yrrphi*rA2*phi  +  ... 

Yrphiphi*r*phiA2  -  (m  +  mx)*u*r; 


K  =  -  zR*(l-t)*T*cos(delta_pod)  +  Kv*v  +  Kr*r  +  Kp*p  +  Kphi*phi  + 
Kvvv*vA3  +  Krrr*rA3  +  ... 

Kvvr*vA2*r  +  Kvrr*v*rA2  +  Kvvphi*vA2*phi  +  Kvphiphi*v*phiA2  + 
Krrphi*rA2*phi  +  ... 

Krphiphi*r*phiA2  +  mx*lx*u*r  -  W*GM*phi; 


N  =  0.5*(l-t)*T*cos(delta_pod)  +  Nv*v  +  Nr*r  +  Np*p  +  Nphi*phi  + 
Nvvv*vA3  +  Nrrr*rA3  +  ... 

Nvvr*vA2*r  +  Nvrr*v*rA2  +  Nvvphi*vA2*phi  +  Nvphiphi*v*phiA2  + 
Nrrphi*rA2*phi  +  ... 

Nrphiphi*r*phiA2; 


%  Dimensional  state  derivatives 

xdot  =[  X*(VA2/L)/m  I  I 

-(-m33*m44*Y+m32*m44*K+m42*m33*N)/(m22*m33*m44-m32A2*m44- 

m42A2*m33)*(VA2/L) 

(-m32*m44*Y+K*m22*m44-K*m42A2+m32*m42*N)/(m22*m33*m44- 
m3  2  A2  *  m44  -m42  A2  *  m3  3  )  *  ( VA2/L  A2) 

(-m42*m33*Y+m32*m42*K+N*m22*m33-N*m32A2)/(m22*m33*m44- 
m3  2  A2  *  m44  -m42  A2  *  m3  3  )  *  ( VA2/L  A2) 
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p*(V/L) 

cos(phi)*r*(V/L) 

(cos(psi)  *u-sin(psi)  *  cos(phi)  *  v)*  V 
(sin(psi)*u+cos(psi)*cos(phi)*v)*V 
deltadot 
ndot  ]; 

D.  V07SIM 

%  SIM.m 

% 

%  Main  program  for  contship  simulation 

% 

%  NTH  1994,  Trygve  Lauvdal 
clear 

v07_data  %  set  model  data 

load  v07_data  %  load  model  data 

N  =  15000;  %  number  of  samples 

table  =  zeros(N,length(x)+length(u));%  storage  vector 
disp(’Simulating’); 

for  i  =  1:N+1, 
table(i,:)  =  [  x' ,  u']  ; 
u  =  v07_cont(i*h,x,u); 
xdot  =  v07_rudder(x,u’); 
x  =  x  +  h*xdot ; 
end 

E.  V07POD  PLOT 

function  plotdata  =  V07pod_plot(  table, h,N) 
elf 


%  START  OF  MAIN  LOOP 
%  data  storage 
%  input 

%  xdot(i)  =  f(x(i),u(i)) 

%  euler  integration  x(i+l) 

%  END  OF  MAIN  LOOP 


figure  (1) 
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x  =  table(:,7); 
y  =  table(:,8); 
plot(x,y,'r') 

grid  on,  xlabel(’X  [m]’),  ylabel(’Y  [m]’) 
title(’Position  in  the  xy  Plane’) 
hold  on 

figure(2) 
time  =  (0:N)*h; 
u  =  table(:,l)*  1.94384449244; 
v  =  table(:,2)*  1 .94384449244; 

U  =  sqrt(u.*u  +  v.*v); 
plot(time,U,’r’), 

grid  on,  xlabel(’Time  [s]’),  ylabel(’U  [knots]’) 
title(’Total  speed  U’) 
hold  on 

figure(3) 

plot(time,v,’r'), 

grid  on,  xlabel(’Time  [s]’),  ylabel(’v  [knots]’) 
title(’Sway  Velocity  v’) 
hold  on 

figure(4) 

phi  =  table(:,5)*(  180/pi); 
plot(  time, phi, ’r’) 

grid  on,  xlabel(’Time  [s]’),  ylabel(’\phi  [deg]’) 
title(’Roll  Angle  \phi’) 
hold  on 
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